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We study the properties and behaviour of the quasi-pseudospherical and quasi-planar Szekeres models, 
obtain the regularity conditions, and analyse their consequences. The quantities associated with "radius" 
and "mass" in the quasi-spherical case must be understood in a different way for these cases. The models 
with pseudospherical foliation can have spatial maxima and minima, but no origins. The "mass" and 
"radius" functions may be one increasing and one decreasing without causing shell crossings. This case 
most naturally describes a snake-like, variable density void in a more gently varying inhomogeneous 
background, although regions that develop an overdensity are also possible. The Szekeres models with 
plane foliation can have neither spatial extrema nor origins, cannot be spatially flat, and they cannot 
have more inhomogeneity than the corresponding Ellis model, but a planar surface can be the boundary 
between regions of spherical and pseudospherical foliation. 



I. INTRODUCTION 

The Szekeres metric is important because, as a model 
with 5 arbitrary functions, it exhibits features of non-linear 
gravitation that less general models cannot. It is an ex- 
act inhomogeneous solution of the Einstein field equations 
(EFEs) that has a realistic equation of state (dust) suitable 
for the post recombination universe, it has no Killing vec- 
tors. It is necessary to pay more attention to models with 
little symmetry in order to better understand all the fea- 
tures and possibilities of General Relativity, and therefore to 
better model the structures of our universe. 

Although there have been a number of papers that in- 
vestigate the Szekeres metric generally 0, 0, [H, @|, and 
several papers that investigate the quasi-spherical case in 
particu lar d @, 0, H, H, Hfl EH, H2] ■ there have been none 
that specifically look at the quasi-pseudospherical and quasi- 
planar cases. This is probably because we have a good un- 
derstanding of spherical gravity from Newtonian theory, and 
so relativistic analyses of spherically symmetric metrics were 
easily developed. Without spherical symmetry, or a slight 
variation of it, familiar relationships, such as that between 
the mass inside a sphere and the gravitational potential, 
do not apply, so it is much more difficult to interpret the 
equations physically. 

We here set out to improve our understanding of the 
quasi-planar and quasi-pseudospherical models, and thus 
enhance their usability, by analysing their physical and ge- 
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ometric properties. The main challenge is to develop a re- 
interpretation of quantities such as "radius" and "mass" 
that cannot retain the meaning they have in (nearly) spher- 
ical models. 

The appearance of the first paper [l3[ to produce an ex- 
plicit model using the quasi-spherical Szekeres metric, that 
of a void adjacent to a cluster, and to plot the evolution 
of its density, is an encouraging development. If the other 
Szekeres cases are sufficiently well understood, explicit mod- 
els can be produced from these too. 

Our methods below are to (a) analyse how the metric 
functions affect the geometry, the matter distribution, and 
the evolution, (b) derive regularity conditions on the metric 
for well behaved matter, curvature and evolution, (c) com- 
pare with other metrics that have planar and pseudospheri- 
cal symmetry, and (d) produce one or two simple examples. 



II. THE SZEKERES METRIC 



In this section we will present the metric and its basic 
relationships, but we will refrain from any physical inter- 
pretation, reserving that for a later section. Once all the 
features and properties of the model are established, we 
will collect the results, discuss the meaning of the various 
functions, and attempt an interpretation of the model. 

Our notation is that of [14|, for which this is a follow up. 
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The LT-type Szekeres metric [H, [ill 1 is: 



or elliptic, / < 0: 



ds 2 



-<2f 



(e + /) 



£ 2 



(2.1) 



where ' = d/dr, e = ±1,0 and / = f(r) > —e is an 
arbitrary function of r. The function E is given by 



e = 0,±l , 



g-Q 
5 



(2.2) 



where 5 = S(r), P = P(r), and Q = Q(r) are arbitrary 
functions. In the original parametrisation of Szekeres, E 
had the form 

E{r,p, q) = A( P 2 + q 2 ) + 2B lP + 2B 2 q + C , (2.3) 
where 2 

A 1 B -- P B - Q 
A -2S> Bl -^S' B2 ~^S> 

p2 , f)2 , c2 

C= 25 ' 4(AC-^-5l) = e. (2.4) 

The function R = R(t,r) satisfies the Friedmann equation 
for dust 



(2.5) 



where ' = d/dt and M = M(r) is another arbitrary function 
of coordinate r. It follows that the acceleration of R is 



R 



-M 



R 2 



(2.6) 



Solving (|2.5p , the evolution of i? depends on the value of 
/; it can be: 
hyperbolic, / > 0: 



M 

R = — (cosh?? — 1) 



(sinhr? — rj) 
parabolic, / = 0: 



/ 



f 3/2 <j(t-a) 
M 



R= / 9M(f-o) 



2\ 1/3 



(2.7) 



(2.8) 



(2.9) 



1 In Ref. [Tal this family of the Szekeres solutions is called the 
P' + family. 

2 In the original parametrisation of Szekeres, the e is an arbitrary 
function of r. If nonzero, this function can be scaled to +1 or —1 
by the rescalings of the other functions: R = yfejii, E = y/\e\E, 
f = |e|/. The scalings cannot change the signs of e and of /. 



R 



(ji — sin rj) = 



M 

FT) 



(1 — cosry) , 



{~fY' 2 o{t-a) 



M 



(2.10) 



(2.11) 



where a = a(r) is the last arbitrary function, giving the local 
time of the big bang or crunch R = and a = ±1 permits 
time reversal. More correctly, the three types of evolution 
hold for / jM 2 ! 3 >, =, < 0, since / = at a spherical type 
origin for all 3 evolution types. The behaviour of R(t,r) is 
identical to that in the Lemaitre-Tolman (LT) model, and 
is unaffected by (p, q) variations. 

The 6 arbitrary functions /, M, a, S, P and Q give us 
5 functions to control the physical inhomogeneity, plus a 
choice of the coordinate r. Note, however, that in the case 
e = we are free to redefine the functions R, S, f and M 
as follows: 



(R, S, f, M) = {aR, S/a, a 2 /, a 3 M), 



(2.12) 



where a = a(r) is an arbitrary function, and the form of 
the metric, the density and the evolution equations will not 
change. In particular, we can choose a so that 5 = 1. 

The density and Kretschmann scalar are functions of all 
four coordinates 



8np = G tt = 



2(M' - ZME'/E) 
R 2 (R' - RE'/E) ' 



(2.13) 



K = R a0lS R aPl & = (8tt) 2 



- 3 0P 



-0PP 



Zp 2 



where 



6M 

8. p = w 



(2.14) 



(2.15) 



is some kind of mean density. For all p and 0p we have K. > 
0, but assumptions of positive mass and density require p > 
and 0p > 0. The flow properties of the comoving matter 
were given for any e value in [14j. For further discussion of 
this metric see [15l, Il6|. 

In the following, we will call the comoving surfaces of 
constant r "shells", and paths that follow constant p &i q 
will be termed "radial". We will use the term "hyperbolic" 
to describe the time evolution for / > 0, and "pseudospher- 
ical" or "hyperboloidal" to describe the shape of the (p, q) 
2-surfaces when e = — 1. To make it clear the shells are 
quite different from spheres, we will call r the "p-radius" 
or "h-radius", R the "areal p-radius" or "areal h-radius', 
and M the "p-mass" or "h-mass", in the planar or pseudo- 
spherical cases, respectively. However, we will use "radius" 
generically when more than one e value is considered. 
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A. Singularities 

The bang or crunch occur when t = a or t = 
2?rA//(-/) 3 / 2 + a, which makes R = and both p and JC 
divergent. Shell crossings happen when surfaces ("shells") 
of different r values intersect, i.e. R' = RE 1 / E and 
M' ^ 3ME'/E. Also p but not K. passes through zero 
where E'/E exceeds M'/3M. 



B. Special Cases and Limits 

The LemaTtre-Tolman model is the spherically symmetric 
special case e = +1, E 1 = 0. 

The Ellis metrics result as the special case E' = 0; 
they are the LT model and its counterparts with plane and 
pseudospherical symmetry. 

The vacuum case is (M' — 3ME' /E) = 0, which implies 
E' = M' = = S' = P' = Q'. For M ^ this gives 
pseudospherical and planar equivalents of the Schwarzschild 
metric [18] (see section fVI Bj) . 

The null limit is obtained by taking / — * oo after a suit- 
able tranformation. In this limit the 'dust' particles move 
at light speed [HHJ 

and the metric becomes a pure radi- 
ation Robinson-Trautman metric of Petrov type D (see [21] 
eq (28.71) with (28.73)). 

The Kantowski-Sachs (KS) type Szekeres metric is in fact 
a regular limit of the LT type Szekeres metric fl5l l20j|. 



C. Basic Physical Restrictions 

1. In order to keep the metric signature Lorentzian we 
must have 



e + f > , 



and in particular 



e + />0 and R'-^jLQ, 



E 



while 



e + / = where R' 



RE' 
~~E~ 



(2.16) 



(2.17) 



(2.18) 



Clearly, pseudospherical foliations, e = — 1, require 
/ > 1, and so are only possible for regions with hy- 
perbolic evolution, / > 0. Similarly, planar foliations, 
e = 0, are only possible for regions with parabolic or 
hyperbolic evolution, / > 0; whereas spherical folia- 
tions are possible for all / > — 1. 

2. We require the metric to be non degenerate & non 
singular, except at the bang or crunch. For a well 
behaved r coordinate then, we need to specify 



oo > 



(R' - RE'/Ef 
(e + f) 



> 



(2.19) 



Whilst failure to satisfy this may only be due to bad 
coordinates, there should exist a choice of r coordi- 
nate for which it holds. 

3. The density must be positive, and the Kretschmann 
scalar must be finite, i.e. 

M 1 - ZME'/E „ 
°° > Ri-RE</E (Z20) 

4. We assume 

R > , M > and S > . (2.21) 

The sign of S, and hence of E can be flipped without 
changing the metric, but S = is not acceptable. 

5. The various arbitrary functions should have sufficient 
continuity — C 1 and piecewise C 3 — except possibly 
at a spherical origin. 



D. 3-spaces of constant t 

It is known from [22| that when e = +1 these 3-spaces 
are conformally flat, and it is easy to verify, using Maple [24] 
and GRTensor [25| that the Cotton- York tensor is zero for 
all e (see [l5|], section 19.11, exercise 19.14, and theorem 
7.1). 

Calculating the Riemann tensor for the constant t spatial 
sections of (12.1ft, we find 
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n-pqpq — 
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(2.22) 
(2.23) 

(2.24) 



ijkl 



(2.25) 

where equations (2)-(5) of [2^| have been used, and the 
other curvature invariants are linearly dependent on these. 
The flatness condition 3 R a bcd — requires 



e^0 
e = 



/ = 
R' = E' 



/' = / = 



(2.26) 
(2.27) 



and the latter is only possible as a limit, or as a Kantowski- 
Sachs type Szekeres model [23]. Interestingly, (e + /) does 
not enter any curvature invariants, and they are all well 
behaved if / = 0. The 2-spaces of constant t and r have 
Ricci scalar 



2 R _ 2e 



(2.28) 
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E. General properties of E(r,p, q) 



and if x — ► ±00 



From (|2.2p we see £7 has circular symmetry about p = P, 
q = Q, which is a different point for each r. The E = 
locus 



(p - P) 2 + (g - Q) 2 - -eS 2 , 



(2.29) 



only exists if e < 0, and is clearly a circle in the p-q plane, 
with E > on the outside, but becomes a point p = P, 
q = Q if € = 0. We have 



E 1 



5 



p'- 



(2.30) 



so the E' = locus is also a circle in the p-g plane, since it 
can be written 



p-P , P'\ 2 , ^ g -Q , Q'V (P') 2 + (Q') 2 



5 



5" 



+ e . 



(2.31) 



With e > 0, this locus always exists, and with e = — 1 it 
only exists if 



y\2 



(s'Y < {p'Y + (q') 



(2.32) 



with the radius of this circle shrinking to zero as the equality 
is approached. Since, if they exist, the distance between the 
centres of these two circles never exceeds the sum of their 
radii 

S_ 



V(P') 2 + (Q') 2 < 

(|5'|V=i +#FW+W) 

(2.33) 

they always intersect, and the intersection points are 



p-P _ eP'S' ± QW-t{( p ') 2 + {Q') 2 +<S')' 2 } 



S 



(P') 2 + (go 2 



g-Q _ eQ'S' T PV-e {(PQ 2 + (Q 1 ) 2 + e(S') 2 } 



(P') 2 + (QO 



l\2 



(2.34) 



To see how E 1 / E affects the metric and the density, we 
write x = E'/E. Then in the metric (|2.ip , g rr is a decreas- 
ing function of x provided x > R' /R, while for the density 
(|2.13|> we have 



87TP 



6M (M'/(3M) - x) 
(P'/P - x) ' 



(2.35) 



8?rp 



6M 

1?3" 



(2.37) 



Therefore at given 7- and i values, the density varies mono- 
tonically with x = E'/E, but the sign of the numerator may 
possibly change as R evolves. If x can diverge, p approaches 
a finite, positive limit. 
The metric component 



{dp 2 + dq 2 
E 2 



(2.38) 



is a 2-d surface of constant unit curvature, that is a pseu- 
dosphere 3 , a plane, or a sphere in Riemann or stereographic 
projection: 

e = -1 , E > : ^~ ^ = coth f cos(0) , 

(i^=coth(Q S in(0), (2.39) 
e = -1 , P < : g = tanh (|) cos(^) , 



^Ql = tanh (0 sin(0) , 

e = : ^ = ff] cos ( 0) 



(<?-<?) 



sin(0) , 



(p-P) 

e = +1 : either — — = cot I - ) cos(</>) 

— s — = UJ ^ 



(2.40) 



(2.41) 



(2.42) 



or 



tan ( - ) cos(<£) , 



= tan (J), into). 



(2.43) 



The projections are illustrated in figs[TJ[3] and the 9-to-p 
transformations (at = 0) are shown in fig [U In these 
diagrams, the parametric equations for spheres and right 
hyperboloids are 

x = a sin cos (f> , y = a sin 9 sin <p , z — a cos , 

< < 7T , < (/> < 2tt , (2.44) 

x = a cos <f> , y = a sin , z = a0 , 

< < 00 , O<0<2tt, (2.45) 

x = a sinh cos , y = a sinh sin <^> , z = a cosh , 
- 00 < < 00 , < </> < 2?r , (2.46) 



so that 



3 The hyperbolic equivalent of a sphere is a right hyperboloid of 
revolution, often called a pseudosphere. 
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FIG. 1: The Ricmann projection from (9,0) to (p, q) coordi- 
nates for pseudospheres (e = —1). The projections of the two 
sheets require different formulae, one is shown as solid grey 
lines, the other as dark dashed lines. The 45° asymptotes 
that divide the projections of the two hyperboloid sheets are 
shown as dot-dash lines. This and the next 2 diagrams show 
only the <j> = 0,7r section, i.e. the q = Q section. For the full 
projection, they should be rotated around the z axis and the 
q dimension added. 



where the former gives the entire sphere minus one point, 
but the latter gives only one sheet of the hyperboloid 4 . No- 
tice that, with 6 &l (f) ranging over the whole sphere, each 
of the spherical transformations f|2 . 39fl & (j2.40[) covers the 
entire p-q plane. (In fig [3] only the range < 9 < ir/2, 
<f> = 0,7r has been shown for each.) In contrast, BOTH of 
the pseudospherical transformations (j2.42fl &l (j2.43|) . with 
< 9 < oo, are required to cover the entire p-q plane once, 
each transformation mapping one of the hyperboloid sheets 
to the p-q plane. To distinguish the sheets, we choose 9 to 
be negative on one and positive on the other. In the planar 
case, the Riemann projection can be considered an inversion 
of the plane in a circle, which is hard to illustrate, or as in 
fig [2] a mapping of a semi-infinite cylinder to a plane. 

One might suspect that the two regions of the (p, q) plane 
on either side of E = simply provide a double covering of 
the same surface, but this is not the case. For the double- 
sheeted hyperboloid at a single r value, the two sheets 
are isometric to each other, the isometry transformation 

(p,q) -> (pW) is 



P = Po 



S 2 (p'-P Q ) 



(p> - P ) 2 + (q> - Q y 



4 Thus, with e = — 1, each constant r "shell" seems to be a hy- 
perboloid with two "sheets" . It will be determined later whether 
both these sheets arc needed or even allowed. 



pv \ / 

i. \ 
i \ » 
i \ ^ 
i s ' 

-P) = S(2/9)j\ \ \ 

1 s \ \ i 




(p - P) = S (2/6) 


: 


l/J 




> 






X 


1 s / / T 

< / / / 

i' / / 

1 ' ; 
1 ' ; 

y i 

■ i 

i / 







FIG. 2: The Riemann projection from (8, <f>) to (p, q) coor- 
dinates for semi-infinite cylinders (e = 0). Section IV CI gives 
the projection as an inversion of the plane in a circle, which 
can't be illustrated as above. The diagram here shows the 
projection of a cylinder, with 8 increasing along the length 
of the cylinder, and each half cylinder maps to the full (p, q) 
plane (with the same formula). 




S COt(6/2) 



FIG. 3: The Riemann projection from (8, <t>) to (p, q) coordi- 
nates for spheres (e = +1). Each of the two possible projection 
formulae maps the full sphere to the plane, but only half of 
each is shown, one as solid grey lines, the other as dark dashed 
lines. 



Qo 



S 2 (q'~Qo) 



[p> - P ) 2 + (q> - Q ) ; 



(2.47) 



where (Sq, Po,Qo) are the values of P, Q and S in that 
hyperboloid. However, for a family of hyperboloids im- 
mersed in a Szekeres spacetime, the transformation (|2.47f) 
will change the values of the functions (P, Q, S) in all other 
hyperboloids, and will not be an isometry. Thus, the two 
sheets are distinct surfaces in spacetime. 
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Thus the factor e determines whether the p-q 2-surfaces 
are pseudospherical (e = —1), planar (e = 0), or spherical 
(e = +1). In other words, it determines the shape of the 
constant r 2-surfaces that foliate the 3-d spatial sections of 
constant t. The function E determines how the coordinates 
(p,q) map onto the 2-d unit pseudosphere, plane or sphere 
at each value of r. Each 2-surface is multiplied by factor 
R = R(t,r) that is different for each r and evolves with 
time. Thus the r-p-q 3-surfaces are constructed out of a 
sequence of 2-dimensional spheres, pseudospheres, or planes 
that are not arranged symmetrically. Obviously, for e < 
the area of the (t = const, r =const) surfaces could be 
infinite, but in the e = +1 case it is 4irR 2 . 



FIG. 4: The relation bewtween p & 9 for each of the Szekeres 
types. Only (p — P) > and 9 > is shown, as rotating <j> 
completes the mapping. The dark dot-dash line is for e = +1, 
the pale solid line is for e = 0, and medium dashed line is for 
e = -1. 
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FIG. 5: An inversion with respect to the circle IC centred at 
CI maps an infinite straight line into a circle that passes through 
CI. The straight lines l\, I2 and £3 are mapped into the circles Ci , 
C2 and C3, respectively. The straight line I4 that passes through 
CI is mapped onto itself, i.e. the image-circle has then an infinite 
radius. A strip between two parallel straight lines is mapped into 
the crescent-shaped ring between their image-circles. The ring has 
finite surface area except when one edge of the strip passes through 
CI. 



It is a property of the Riemann projection that for e > 0, 
circles in (p,q) map to circles in (9,<fi). Constant </> lines in 
(6, <f>) (that obviously pass through 8 = 0) map to straight 
lines through p = P, q = Q. Circles in (8,<fr) that pass 
through 8 = 0, map to straight lines in (p,q). See fig. [5] 
for an example with e = 0. For e = the projection is just 
an inversion of the plane in the circle of radius V25- 



III. THE EFFECT OF e AND E 

We here analyse the role E plays in these models, and 
contrast it with the e = +1 case, in which E 1 /E creates 
a dipole variation around the constant (i, r) 2-spheres. We 
omit some of the detail below because very similar calcula- 
tions were done in [14]. We assume S > 0. 



A. Pseudospherical foliations, e = — 1 



Transforming (j2.2j) and its derivatives using (|2.39p and 
(|2.40p , and putting e = — 1 , we get 



E 



vS 



cosh 8 — v 



(3.1) 



, S' cosh 8 + sinh 6»(P' cos <t> + Q' sin cj>) 

-k = r~2i > l 3 - 2 ) 

cosh 8 — v 



E" 



S" cosh 8 + sinh 8{P" cos + Q" sin < 



(cosh# — v) 



S' cosh 8 + sinh 8(P' cos <f> + Q' sin < 



(cosh# — v) 



{{S'f - {P'f {Q'f) 
S 



(3.3) 



where v — +1 when E > 0, when E = and —1 when 
E < 0. The E = circle corresponds to 6 — * ±00, and 
its neighbourhood represents the asymptotic regions of the 
two sheets. It is clear that curves and regions that intersect 
the E = circle must have infinite length or area, since 



L = J ^\/dp 2 + dq 2 ds , A 



R 2 

-^dpdq (3.4) 



The locus E' = for all E is 

S' cosh 8 + P' sinh 8 cos (j> + Q' sinh 9 sin cj> = . (3.5) 
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Writing z = cosh#, y — sinh6*sin0, x = sinh^cosi/) as 
the parametric locus of a unit right hyperboloid centered 
on (0,0,0) in flat 3-d space, we find (|3.5p becomes S'z + 
P'x + Q'y = which is a plane through (0,0,0), so E' = 
is the intersection of a plane with a right hyperboloid. In 
fact, (|3.5[) is a geodesic of the p-q 2-space, as shown in 
appendix [A] 

We can write the E' = locus as 



tanh 9 



-S' 



-d 



P' cos tp + Q' sin cf> cos( 



(3.6) 



where 



S" 



V(P') 2 + (Q') 2 ' 



tan(0 o ) = ^ (3.7) 



so obviously a solution only exists if ()2.32jl holds, and only 
for 



| cos(0 - O )| > d . 
From (I3T2D and (|3~T1> we find 



(3.8) 



5" cosh + sinh0(P' cos </> + Q' sin 0) 
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thus £"/_E = constant implies S'z + P'x + Q'y — Sx con- 
stant, which is a plane parallel to the E' = plane. The 
location of the extrema of E'/E are found as follows 

d{E'/E) _ smh g(P' sin </>-Q' cos 0) 



P' 



v/(P') 2 + {Q'f 



(3.11) 



d{E'/E) 



de 



^ S' sinh 6> + cosh 6(P' cos (/> + Q' sin <j>) 



tanh^p = — 



S 

P' cos 4> e + Q' sin ( 
J' 



5" 



cosh 9 e = e 2 



S' 



vW - (^') 2 - (Q'Y 



(3.13) 



where ei = ±1 and €2 = sign(S"). The extreme value is 
then 



V(S')2 - (PQ 2 - {Q'Y 
e 2 y . (3.14) 



extreme 

and these extrema only exist at finite 9 if 

{S'f > (P'f + {Q'f 



(3.15) 



which is the opposite of (|2 . 32fl : so on a given constant r 
shell, either E' = exists, or the extrema of E'/E exist, 
but not both. Notice that when ()3.15p holds, then E' does 
not change sign on a given sheet, it is fixed by v and the 
sign of S' . It follows from (|3.9p that this extremum is a 
maximim where E'/E is negative, and a minimum where 
E'/E is positive. Thus, for each constant r hyperboloid, 
on the sheet with ES' < (i.e. vt 2 = -1), E'/E has a 
positive minimum and goes to +00 as \0\ — > 00, while on 
the sheet with ES' > 0, E'/E has a negative maximum and 
goes to —00. The maximum and minimum are at opposite 
poles in the sense that {9, (f>) — * {—9, 4> + tt) maps one into 
the other, and indeed it maps E'/E to —E'/E. We now 
specify that 9 < on the E < sheet (see below (12. 46ft ). 

From the foregoing considerations, if {S') 2 > (P') 2 + 
{Q') 2 , then E'/E is the pseudospherical equivalent of a 
dipole, having a negative maximum on one sheet and a pos- 
itive minimum on the other, but diverging in the asymptotic 
regions of each sheet near E = 0. 

We see in the metric (|2.ip that RE' /E is the correction 
to the separation R' , along the r curves, of neighbouring 
constant r shells, meaning that the hyperboloids are cen- 
tered differently and are "non concentric". In particular 
RS'/S is the forward {9 = 0) displacement, and RP'/S & 
RQ' / S are the two sideways displacements {9 = tt/2, cf> — 
0) & {9 = it/2, (p = tt/2). The shortest radial distance 
is where E'/E is maximum. (From a given point (p,q) 
on a given r shell at constant t, the shortest distance to 
an infinitesimally neighbouring r shell must be along an or- 
thogonal curve, i.e. along constant p and q.) 



B. Planar foliations, e = 



Transforming (|2.2p and its derivatives using (|2.4ip and 
putting e = , we get 



F- 2S 



(3.16) 



E , = 2{S' + 9{P' co^ + Q'sm^)) 



E" 



+ 4 



e 2 



2{S" + 9{P" cos c/) + Q" sin0)) 
J 2 



S'\ fS ' + 9{P' cos + Q' sine 

02 



(PT-_(Q') 
S 



'\2 



(3.18) 



The E = locus has shrunk to the point p = P, q = Q, 
but still corresponds to the asymptotic regions of the plane, 
6 = 00. The locus £" = is 



5' + P'9 cos + Q'(9sin0 = 



(3.19) 
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Obviously, (|3.19p is a geodesic of the p-q 2-space. We 
can write the E' = locus as 



-5" 



(3.20) 



P' cos + Q' sin cos(0 — 0o) 
where (|3.7p defines d and cj)Q, and evidently it exists provided 



and (P' ^ or Q' ^ 0) 

From (gUZD and (|3~T6]l we find 

E' _ S' + 0{P' cos cj) + Q' sin 0) 
~E ~ S ' 



(3.21) 



(3.22) 



Thus there are no extrema of E 1 / E, and it extends to both 
±oo, though for fixed 8, = 0o ± 7r gives the line of maxi- 
mum and minimum E'/E. The behaviour found here can- 
not really be termed a dipole. 

As before, RE' /E is the correction to the "radial" sepa- 
ration R' of neighbouring constant r shells, and the above 
indicates that adjacent shells are planes tilted relative to 
each other, with 0o being the direction of maximum tilt, 
but if E'/E = constant they are parallel. 



IV. REGULARITY 

A. Pseudospherical and Planar "Origins", 

For spherical foliations, e = +1, if r = is an origin, 
then R(t,0) = for all t, and such origins are well under- 
stood. The conditions on the arbitrary functions that ensure 
a regular origin were given in [14]. Specifically, the density, 
curvature and evolution of R are all well behaved if 



S-R 71 



M ~ R 3 
P~R n 



f~R 2 
Q ~ R r 



n > 0. 



(4.1) 



We note that the derivation of these conditions does not 
depend on the value of e. Therefore one immediately asks 
whether such a locus is possible for pseudospherical and 
planar foliations. 

Now by (|2.16p we must have / > — e for a Lorentzian 
signature, so for e = — 1 models, / — > is not possible. 
Therefore an "origin" is not allowed for pseudospherical fo- 
liations. 

For planar foliations, e = 0, / — > is not impossible. By 
(|2.19p we expect 



R' 



lim g r 

i — >ro 



lim ■ 

r^ro 



1 



BE' 
R'E 



f 



(4.2) 



to be finite and non-zero, and from (79) and (84) of [14] 
we know RE' /(R'E) is not divergent. So, to keep g rr well 
behaved in this limit, we require R'/yff to be finite and 
non-zero, and by ()4.ip this implies 



R'~^f~R 



R 



b constant, (4.3) 



while the "radial" distance is 

s = / Jg r r dr ~ r 



(4.4) 



In other words, R only asymptotically approaches zero. 
Therefore there is no real origin, but R, M and / can 
asymptotically approach zero. (See figure [9]) 



B. Conditions for No Shell Crossings 



For p to be positive, (j2~13|l shows that (M 1 - 3ME'/E) 
&l (R' — RE'/E) must have the same sign. We now consider 
the case where both are positive. Where (M'-3ME'/E) < 
and (R' — RE' /E) < we reverse the inequalities in all 
the following. 



1. Pseudospherical foliations, e 



The inequality 



(M' - 3ME'/E) > 



(4.5) 



must hold for all possible p &l q, and at every r value. If 
(|2.32p holds so that there is an E' = locus on each hy- 
perboloid sheet, then E'/E varies between ±oo, diverging 
in the asymptotic regions of each sheet, so the density in- 
evitably goes negative in some regions of every constant r 
shell. If however (|3.15p holds, so there are finite extreme 
values for E'/E but no loci where E'/E = 0, then on 
the sheet with < (E'/E) mi „ < (E'/E) < oo, |g3J is 
violated over all of the sheet, except near the minim um 



if M'/(3M) > (E'/E) mi 
but on the sheet with > (E'/E) n 
is always satisfied if 



(P'Y~(Q'Y/s, 

> (E'/E) > -oo, it 



M' 
3M 



> 



y/(S')2 - (P'Y - (Q'Y 

s 



(4.6) 



It is obvious that pT5]) and (j4~6]) ensure ([43] ), but (]4~6j) 
can only hold for one sheet, and on that sheet it appears 
that negative M' f or R' are not excluded. 

Now consider the time evolution of (R'/R - E'/E) > 0. 
Because of the above, we only need consider the negative 
E'/E sheet, and since e = — 1, only hyperbolic evolution, 
with / > 1, is relevant. The argument proceeds almost 
exactly as in [141 ], except that the dipole term —RE'/E 
is everywhere positive, so it tends to relax the conditions. 
Defining cf>4 = sinh^sinlirj — ??)/(cosh?7 — l) 2 and 05 = 
sinh?//(cosh77 — l) 2 we have 



R ' M 'n a>\+? f 3 ^ 

I=m (1 ^ 4) + 7^ 4 



1 1 f 3/2a \ 



(4.7) 



Because (1 — 04), (304 — 2) and 05 are always positive, but 
evolve differently with ?/, this argument shows that to avoid 
shell crossings we require 



a' < 



(4.8) 
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and 



2/ E 



>0 , 



(4.9) 



and the latter takes its strongest form at the maximum of 

E' jE, so 



f V(S') 2 - (P') 2 - {Q'Y 

2/ " 



S 



(4.10) 



To confirm (pTT5|> . ([4~6]l . ([4~8]l and ([4~TUll are sufficient, 
we use equations (99) and (100) of [14J, write X = 
\(E'/E) max \ so gSJ and pTOj) become M'/(3M) = 
— X + a and /'/ (2/) = — X + with a and /3 non-negative, 
and thus obtain (|4.7p again and 



i?' f 3 / 2 a' 
- = -X + 3a(l - 4 ) + /3(30 4 - 2) - ^^-05 



as required. This also means 

R' M' fo . tnl n ^ f/ 2 a' 
> (0 - a)(36 4 - 2) - 5 



(4.11) 
(4.12) 



(4.13) 



so the numerator of (|2.36|) is negative for all r\ if there are 
no shell crossings and 



/' M' 
2f ~ 3M 



(4.14) 



otherwise it can change sign from negative to positive as r\ 
increases or if a' = it goes from zero to negative. 

Thus we see that only one of the hyperboloid sheets can 
be free of shell crossings, and it must have a minimum in 

(R' - RE' )E) and (M' - ME'/E). 



2. Planar foliations, e = 



By (|3.22p and the discussion following (|3.19p we have 

R! E' R' S' + 9{P' cos cj) + Q' sine 



R E R 



S 



, (4.15) 



so adjacent shells are like tilted planes and inevitably they 
must intersect on the straight line 



-(S' + SR'/R) 

o = 

P' cos 4> + Q' sin q 

creating shell crossings, except when 

P' = , Q' = , 

and 

R!_ -S' 
R - S ' 



(4.16) 



(4.17) 



(4.18) 



Condition (|4.17p ensures the shells are parallel, while (|4.18p 
can be converted to 



R' n 
R^° 



(4.19) 



because S can be absorbed into other functions, as shown 
in equation (|2.12p . Effectively then we require 



S' = P' = Q' = E' = 



(4.20) 



and the remaining conditions follow exactly as in |14j or 

The no shell crossing conditions for e < are summarised 
in Table 1. It is a continuation of Table 1 in Sec. VI of Ref. 
UM . which summarised those conditions for e = +1. 



C. Regular Maxima and Minima 

We already know that spherical foliations can have reg- 
ular extrema r = r m , where R'(t,r m ) = 0, and we con- 
sider this possibility for other e values. The case of both 
(M' - ZME'/E) & {R 1 - RE' jE) being zero may occur 
momentarily at isolated locations as R evolves, but for a 
given r = r m , the no shell crossing considerations give 

R' = M' = f = a' = S' = P' = Q' = , (4.21) 

since they must hold at all times, and for all p &l q. In order 
for the metric and the density to have well behaved limits 
as r m is approached, we require 



j^f J^E' IE 

g rr = [J — — > L , < L < oo (4.22) 



AnpR 



2 _ M' - UIE'/E 
R' - RE' )E 



N 



< iV < oo 



(4.23) 



and we obtain all the results of section VI of [14], which 
was done for general e. As noted there, we must replace M' 
with lim r ->r m M' /y/e + f , and similarly for all 6 arbitrary 
functions, in all the no shell crossing conditions; and to 
ensure these limits exist as well as avoid a surface layer at 
r m we require 



f = -e 



(4.24) 



With pseudospherical foliations this just means / = +1 at 
an extremum, but with planar foliations, we already saw in 
section [TV Al that / = is not possible, and can only be 
approached asymptotically, so spatial extrema of R cannot 
occur when e = 0. 
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Table 1. Summary of the conditions for no shell crossings or surface layers. 



e 


R' 


/ 


S' 


M' , /' , a' , P' , Q' 


= -1 


> 


> 1 


ES' > 


(S') 2 > (P') 2 + {Q') 2 

M' s y/{S') 2 -{P') 2 -{Q') 2 
3M — S 

f \/(S') 2 -(P') 2 -(Q') 2 
2/ - S 

a' < 




= 


= 1 


S' = 


M' = , /' = , a' = , 
P' = , Q' = 




< 


> 1 


ES' < 


(5') 2 > (P') 2 + (<2') 2 

M' ^ , ^(S') 2 -(P') 2 -{Q') 2 

3M — ' S 

f / , y/(S') 2 ~(P') 2 -(Q') 2 

2/ - + S 

a' > 


= 


> 


> o 


= 


M' > , /' > , a' < , 
P' = , Q' = 




= 


= 


= 


M' = , /' = , a' = , 
P' = 0, Q' = 




< 


> o 


= 


M ' < , /' < , a' > , 
P' = , Q' = 



D. Density: Extrema, Asymptotics, and Evolution 



Considering the density f|2 . 13[) with hyperbolic and 
parabolic evolution (f27D - ([Z9]) and R'/R given by (|40T|) 
and section V.B.I of [14], we have 



070 

a' = 



R 
R!_ 
R 



oo , 

M' 
3M ' 



P 
P 



M' - ZME'/E 

4ttR 2 R' 
3M 



(4.25) 



AttR 3 



PUT early (4.26) 



at early times, 77 — > 0, while at late times, 77 — > 00, 

W - 3ME'/E 



RJ_ 
R 



L 

2/ 



47ri?3(/V(2/) - 



(4.27) 



Therefore, the effect of E' /E only disappears near a simul- 
taneous bang. 

We saw in section fTTTAl that when e = —1, E'/E acts 
like the pseudospherical equivalent of a dipole, and section 
I IV Bl showed E'/E must be negative but rise to a maxi- 
mum somewhere. By (|2.36p and (|4. 13|) the density p de- 
creases monotonically with E'/E if (|4.14p holds, otherwise 
it can change to a monotonic increase as time passes. So we 
conclude that, on shells where (|4.14p holds, the density is 
minimum where E'/E is maximum and the shell separation 



minimum, and vice-versa 5 . On shells where it doesn't hold, 
the initial density minimum can evolve into a maximum. 
This holds for hyperbolic evolution with any e value. 



V. THE CASE OF e = 

Ironically, the e = case is the most tricky to understand. 
Below we consider the quasi-planar case in two ways; as a 
complete manifold with planar foliation, and as a boundary 
surface between a region having a spherical foliation and 
one having a pseudospherical foliation. 



A. The Quasi-Planar Manifold 

It is difficult to interpret the Szekeres spacetime in which 
all the (p,q) subspaces are flat, even in the limit E' = 0, 
when the spacetime becomes plane symmetric. As seen 
from section [TT"Dl the value / = is not admissible, as it 
makes both the metric and the curvature singular. Thus, 



5 It is amusing to note that the M' = 0, M > 0, case could be 
called a "bare dipole". Of course, this case suffers from shell 
crossings and negative densities. 
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the quasi-planar case does not admit flat 3-dimensional sub- 
spaces, so this case cannot provide a foliation of 3-d Eu- 
clidean space, such as the construction of section IVTT1 



B. The Quasi-Planar Szekeres Metric as a Limit 

We here show that the planar metric can be viewed as 
the limit of the other two at large R. We pay particular 
attention to the limit of the spherical case, with which we 
are more familiar. In spherical coordinates, if R is large (at 
a fixed, finite t), the region near 9 = looks like cylindrical 
coordinates, and in the limit as R diverges, the constant 
r surfaces are effectively planar. In this limit pseudospher- 
ical coordinates also look cylindrical. We need to find a 
transformation that will allow this limit but keep all physi- 
cal quantities well behaved. Let w be a large quantity that 
goes to oo in the limit, then the transformation 



M -> uj a 0M , 



S 
t - 
6- 

V - 

P ' 



0t , r - 



0T) 
0P 



R 



/ -> uj 2 0f , 0a , 

P^0P , Q -» 0Q , 

0r , p — > 0p , q—>0q 

<f> — ► 0(j> , 

-> UJ0R , E — ► lo0E , 



(5.1) 



results in 



E =- 



'p-P 
. S 



s 



0E = 



0S\(0p-0P\ 2 (0q-0Q x2 



0S 



0S 



— ^- = cot f - j cos(0) -> = - cos 



(q-Q) 



cot ( - ) sm(4>) 



1-0Q) 2 . 
^5— = ^ Bm( ^ ) ' 




_ -{S"cos6» + sin6'(P'cos0 + Q'sin ( /))} 
£" -{5' + 6»(P'cos</> + Q'sin^)} 



e + f 
R 2 sin 2 6» d<i 2 -► 0R 



0/ 

2 w/i2 ±2 



(5.5) 

. (5-6) 
(5.7) 



8ir0p 



2 (0M' - 30M0E'/0E) 
0R 2 (0R 1 - 0R0E'/0E) 



(5.8) 



R- 



2M 



20M 

0i? 



+ 0/ , 



(5.9) 



R = ^-(coshr? — 1) — ► 0i? = ^^-(cosb.077 — 1) 
/ 0/ 



t — a = 



M 



(sinhrj — ?/) — > 0t — 0a = 



0Af 

0/ 3 / 2 



(5.10) 

(sinh 077 — 0rf) 
(5.11) 



Thus we have exactly the planar Szekeres metric, with all 
the correct matter content and dynamics. Another way of 
looking at this transformation is that we have effectively 
taken an infinitesimal region near 9 = at finite r and 
blown it up to finite size. Note that / must diverge, so 
while an elliptic model can have infinite R [26], it cannot 
have this limit. 



C. The flat limit of the Riemann projection 



Equations (|2.41[) show that the transformation from the 
(p,q) coordinates to the (0,4>) coordinates involves an in- 
version of the coordinate plane in the circle of radius \/2S: a 
point at a distance u = yj (p — P) 2 + (q — Q) 2 from (P, Q) 
is mapped into a point at a distance 9 = \/2S/u, so that 
the product 9u is the same for all point-image pairs. The 
problem is thus to set up the two other mappings in such 
a way that in the limit of zero curvature (infinite radius) of 
the sphere or hyperboloid they go over into an inversion of 
the plane 6 . 

The characteristic property of the inversion is that the 
inversion circle remains invariant. The first question is thus: 
is it possible that in the other two projections a circle in the 
curved surface is mapped into a circle of the same radius 
in the plane? This must be answered separately for the 
sphere and the hyperboloid, and we now proceed to this 
consideration. 



1. The quasi-spherical model 

For the quasi-spherical model, the projection is illustrated 
in Fig. [5] The radius of the sphere is a, the point on the 
sphere that is being mapped has the polar coordinate 9 
and the projection plane PP is at the distance S from the 
projection pole O. The image-point in the plane is at the 
distance p from the axis. (We set P and Q to zero, since 
their values are unimportant for any one shell.) 



8TTP : 



2 (M' - 3ME'/E) 
R 2 (R> - RE 1 IE) 



6 In terms of the limit of the preceding subsection, we have S — > 
0S/ui and 8 — > 00 /uj as u) — > oo, which does ensure ||2.42| | and 
( f239l go to JZ2TJ. 
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FIG. 6: The Riemann projection of a sphere on a plane. 
If the plane intersects the sphere, the circle of intersection 
has radius b = \J a 2 — (S — a) 2 and it is invariant since it is 
mapped onto itself in the projection. For the projection shown 
here, however, in the limit a — > oo with b held constant, an 
identity mapping results, not an inversion. 




FIG. 7: The Riemann projection of a sphere that goes over into 
the inversion of the plane in the limit a — > oo. The meaning of the 
symbols is the same as in Fig. \E\ The limit is taken in such a way 
that the circle of intersection, of radius b, remains the same as the 
radius of the sphere a goes to infinity. One of the larger spheres is 
shown. If the point in the plane at the distance p from the axis is 
kept constant for all spheres, then its image on the various spheres 
will follow the circle arch pt, and in the limit a — > oo the image 
will land in the plane, at the distance b 2 /p from the axis - i.e. the 
limiting plane undergoes an inversion. This result is derived in the 
text. 



If PP intersects the sphere so that the radius of the inter- 
section circle is b, then points on the sphere left of PP map 
to points on the plane outside that circle, and vice-versa. 
This will become the invariant circle of the inversion in the 
limit. For a given b value, there are two possible locations 
for the plane, S = a ± \Ja 2 — b 2 . The "+" configuration, 
shown in Fig. [His, however, unsuitable for the limit of in- 
finite radius, because the part of the sphere right of PP is 
mapped onto the inside of the circle in the plane, and in the 
limit a-*oowe will not get an inversion, but an identity 
mapping. 

Therefore we proceed to the "— " configuration, 



5' 



sja 2 -b 2 , 



(5.12) 



shown in Fig. [JJ We begin with a sphere of radius a, and 
increase a while moving the center of the sphere to the right 
in such a way that all spheres intersect the plane PP along 
the same circle of radius b. We have tan(#/2) = S/p, and 
so 



2aSp 

D 2 +S 2 



2ap (a — \/ a 2 — b 2 ) 



(a — \J a 2 — 6 2 ) + ; 



(5.13) 



b 2 1 {a + \J a 2 — fo 2 ) and obtain 
2ab 2 p 



(a + y/a 2 - b 2 ) 



( t V 



P ' 



(5.14) 

which is indeed an inversion in the circle of radius b. Fig. 
[7] shows also the trajectory of the projected point on the 
sphere as a — > oo, while p is kept fixed. That trajec- 
tory is a circle of radius (l/2)(6 2 /p — p) with the centre 
at (z,x) — (0, (p + b 2 /p)/2). As should be expected, the 
circle degenerates to a point when p = b, and its radius 
becomes infinite when p — > 0. 

The S of the flat case is actually the b = y/S(2a - S) 
of the spherical case, and b and h correspond to 0S and 09 
in (1511). 



We now apply the identity a 



y/a 2 - b 2 



2. The Quasi-Pseudospherical Model 

As with the spherical case, the hyperboloid and the plane 
must intersect for an invariant circle to exist, and if the 
pole of projection is not placed in the same sheet of the 
hyperboloid as the invariant circle, then in the limit of zero 
curvature an identity instead of an inversion results. The 
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FIG. 8: The Riemann projection of a hyperboloid that goes over 
into the inversion of the plane PP in the limit a — > oo. Three 
hyperbolae (intersections of hyperboloids with the plane of the fig- 
ure) are shown. The parameter a has the smallest value on the 
hyperbola with the leftmost vertex and largest for the hyperbola 
with the rightmost vertex. The calculations in the text are done 
for the middle hyperbola. As a increases, and /3 = a — S increases, 
the hyperbola is shifted right so that the circle of intersection of 
the hyperboloid with the plane PP is always the same and has 
radius b. The curve pt is the trajectory followed by a point on the 
hyperboloids as a — > oo, while the image point in the plane PP is 
kept at the same distance p from the axis. All hyperbolae are right 
(their asymptotes are inclined at 45° to the x-axis); those that look 
wider open are simply magnified. 



case that gives the inversion, is shown in Fig. [5] 

We will increase a — > oo, but will shift the hyperboloids 
so that they intersect the plane of projection PP always 
along the same circle of radius b. Therefore we must have 



S 



\fa 2 



b 2 - 



(5.15) 



(We do not consider S — y/ a 2 + b 2 + a which does not 



lead to inversion.) As expected, in the limit a 
S — > 0. Using this we get 



y/a 2 + b 2 



a 



h 2 - 



h 



oo we get 



(5.16) 



which is an inversion in the circle of radius b. 



D. Joining Spherical and Pseudospherical Foliations at a 
Planar Boundary 

Suppose in an e = +1 Szekeres metric we let the radius of 
the constant r spheres diverge, so they become effectively 
planar at some r value, and similarly in an e = —1 Szekeres 
metric we let the "radius" of the hyperboloids diverge at 
some r value. Then the the two metrics can be joined at 
their planar boundaries, provided we carry out the planar 
limits of (|5.ip - (|5.1ip . as is easily verified by calculating the 
junction conditions. 



VI. COMPARISON WITH ALLIED METRICS 

A. Foliations of the Robertson-Walker Metric 

To better understand the e = — 1 and e = Szekeres 
foliations, we first look at the simplest possible cases — the 
homogeneous ones; i.e the Robertson-Walker (RW) metric 
with planar and pseudospherical foliations. The RW metric 
in standard coordinates is 



ds 2 = -dt 2 +S 2 (t)< 



^_ + r 2 (d$ 2 +sm 2 $d<p 2 )\ , 
1 — kr l 

(6.1) 



and the spherical foliations (e = +1) obtained with f — 
sin(r), f = r and f = sinh(r) for k = +1, k = and 
k = — 1 respectively are familiar. For k = — 1 in particular, 
we have 

r = sinh(r s ) , d = 9 S (6.2) 
-> ds 2 = -dt 2 + S 2 {t){dr 2 s 

+ smh 2 (r s )(d8 2 s + sm 2 (6 s )d^ 2 )} , (6.3) 
-> R s = S sinh(r s ) . fs = + sinh 2 (r s ) , 

M s = M sinh 3 (r 5 ) . (6.4) 

For planar foliations, e = 0, with k = we obtain 

f = \Jr 2 + 9 2 , d = tan" 1 - (6.5) 

r 

-» ds 2 = -dt 2 + S 2 (t){dr 2 P0 + {d6 2 PQ + 9 2 P0 d(f> 2 )} , 

(6.6) 

-> Rpo = S , f P0 = , M P0 = M ; (6.7) 
while with k = — 1 we get 

rp = ln(cosh rs + ei sinhrs cos 8s) , &i = ±1 

(6.8) 

, ei = ±l (6.9) 



ei sinh rs sin 8s 



(cosh rs + 62 sinh rs cos 9s) 
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sinhrg = 



i(e^(^ + l) + e-^) 2 -l 



(6.10) 



ei sine's = 



e 2 cost's = 



(| (e™»(0» - 1) + e-^)) + {9 P e^f 

(6.11) 



l) + e-^)) 



(I (e rp (^ - 1) + e-^)) + {Ope^f 

(6.12) 



-> ds 2 = -dt 2 +S 2 (t){dr 2 P +e 2rp (d8p + 9pd4> 2 )} , 

(6.13) 

-» i?P = 5e rp , / P = +e 2 ^ , M P = M e 3 ^ , 

(6.14) 

For pseudospherical foliations, e = —1 &i k = —1, we find 



r = \/sinh 2 (r#) + cosh 2 (r# ) sinh 2 (#ff ) , 
/ sinh(6lff) \ 



d = tan 



(6.15) 



\tanh(r/f) / 

ds 2 = -dt 2 +S 2 {t){dr 2 H 

+ cosh 2 {r H ){d9 2 H +sinh 2 (9 H )d<j) 2 )} , 

(6.16) 

R H =S cosh(rTf) , f H = + cosh 2 (r H ) , 

M H = M cosh 3 (r H ) . (6.17) 



has R(to,rp) decaying asypmtotically towards zero in one 
direction. 

Each constant r# 2-surface has the geometry of one sheet 
of a two-sheeted right hyperboloid of revolution. The func- 
tion coshrfj has a minimum (at r# = 0), which is sugges- 
tive that a "natural" way to cover such a negatively curved 
manifold with hyperboloids is to have R(to,rn) > 0, but 
going through a minimum. 

In the spherical foliation, M and / are constant on 
spheres, are zero at an origin, and reach a maximum where 
R' = 0. In the pseudospherical foliation, the correspond- 
ing M and / functions are constant on completely different 
surfaces, have a minimum where R' = 0, and are nowhere 
zero. In the planar foliation, M and / are again constant 
on different surfaces, have no origin and no extremum, but 
asymptotically approach zero. Despite the apparently very 
different descriptions, these 3 foliations of the k = — 1 case 
describe exactly the same metric with the same behaviour. 
In each case they obey 



S 2 -^-k 



Anp = 



3Mn 



s 3 



M = 



(6.18) 



where po and Sq are constants. (Due to the homogene- 
ity, the features of R, M &l f such as the origin in (|6.4p 
or the minimum in (|6.17p are not special locations, as a 
transformation could move them to any position.) 

To verify that R cannot go to zero in pseudospherical 
foliations, we write the RW metric in the form 

ds 2 = -dt 2 + S 2 (t){x 2 (r) dr 2 + y 2 (r)(d9 2 + z 2 {9) d(f> 2 )} , 

(6.19) 

and require that it satisfy the EFEs with the usual Fried- 
mann equation (j6. 18|) for S. We find 



Using the Riemann transformations, each of these can be 
converted to Szekeres form, but we are here interested in 
understanding the relationship between different foliations. 
The relationship between (rH,9n), (rp,9p) and (rs,9s) 
for the k = — 1 foliations is illustrated in figO which plots r 
and § as polar coordinates on the plane. (This compression 
of a negatively curved 2-surface onto the plane naturally 
creates some distortion.) 

We now consider time sections t = to of the k = — 1 RW 
model. Each constant r$ 2-surface has the geometry of a 
sphere. The function sinrg has a zero (at r$ — 0), a max- 
imum (at rs = 7r/2), and another zero (at r$ — tt), which 
are features of a closed surface in spherical coordinates. 

Each constant rp 2-surface has the geometry of a plane, 
but in order to embed the plane into a negatively curved 
3-space, it has to bend round so that the circumference 
e rp 8pAcf> does not increase too fast compared with the ra- 
dius e rp 9p. The constant 9p surfaces are horns that flare 
out rapidly, and even bend backwards to stay orthogonal 
to the constant rp planes. The coordinates on each plane 
are magnified by the factor e Tp that is nowhere zero, sug- 
gesting that a plane foliation of a negatively curved space 



y 



\Ja- ky 2 (r) 



A9 2 



+ a z = 



(6.20) 



where a is arbitrary and y(r) is not fixed, except when a — 
and k = 0, in which case 



dy = 

dr 



d9 2 







(6.21) 



and x(r) is not fixed. Clearly, if k = +1 a must be pos- 
itive, so z(9) must be a trig function, but if k = — 1 a 
can have either sign, so z(9) may be a trig function or a 
hyperbolic trig function. Thus y(r) = cosh(r) gives (|6.16p 
above. Notice however that if k = —1 and y(r) goes to 
zero somewhere, such as y(r) = sinh(r), then a cannot be 
negative. 



B. Matching the Szekeres Metrics to Vacuum 

We will now match the general planar and pseudospher- 
ical Szekeres solutions to vacuum metrics. Although the 
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FIG. 9: A section through the spherical, planar and pseu- 
dospherical foliations of the fc = — 1 RW model. The curves 
of constant rs and constant 8s are the thin dotted lines, the 
curves of constant rp are the solid pale arches, the curves of 
constant 6p are the solid pale lines diverging from the left, the 
curves of constant th are the medium dashed vertical lines, 
and the curves of constant 6h are the medium dashed curves 
in the left-right direction. This diagram is distorted because 
a negatively curved 2-surface has been compressed onto a flat 
plane. In fact all three sets of lines are orthogonal. The 3-d 
diagram is obtained by rotating this one about the z axis. At 
the right of the diagram it is evident how the 3 coordinate 
systems are similar near the z axis. This becomes exact as 
z — » oo. 



ify 7 . With e = — 1 we obtain the vacuum pseudospherically 
symmetric metric. 

Note that the vacuum metrics with e < are very differ- 
ent from Schwarzschild's. The Schwarzschild metric is static 
for R > 2m and nonstatic (vacuum Kantowski-Sachs) for 
R < 2m, but the two regions together form one complete 
manifold, as evidenced by the Kruskal-Szekeres extension. 
The two other metrics are globally nonstatic when m > 0, 
as T is a space coordinate and R is time. The Kasner so- 
lution with to > is of Bianchi type I, the pseudospherical 
one is of Bianchi type III (see Appendix [B|) . 

First, we write the vacuum metrics in Szekeres form. Fol- 
lowing the prescription used for the Schwarzschild solution 
(see Exercise 10 in Chap. 14 in Ref. [15]), we transform 
(|6.22p to coordinates defined by observers freely falling in 
the indirection — the plane- and pseudospherically sym- 
metric analogues of the LemaTtre-Novikov coordinates for 
the Schwarzschild solution, see [27| and Sec. 14.12 in Ref. 
(T5j . We substitute 



T = T(t,r), R = R(t,r) 



(6.23) 



and require that in the (t,r) coordinates the component 
g t t of the metric is — 1, while g tr = 0. We solve this set 
of equations for T t and T r , then impose the integrability 
condition T.t r = T >r t- Discarding the trivial case R. r = 0, 
it reduces to 



(6.24) 



where F(r) is an arbitrary function. This is a special case 
of eq. ()2.5p . corresponding to M = to = const and F = f. 
The full solution for T(t,r) is given in [23[. The metric 
(|6.22p in the (t,r) coordinates becomes 



vacuum metrics are very different in each case, the match- 
ing can be solved at one go for all e. Inspired by Bonnor's 
result [6, 7] that the quasi-spherical Szekeres metric in its 
full generality can be matched to the Schwarzschild solu- 
tion, we will verify that the two other Szekeres solutions, in 
their full generality, can be matched to the corresponding 
plane- or pseudospherically symmetric vacuum solutions, re- 
spectively. The planar and pseudospherical analogues of the 
Schwarzschild solution are known, even if not well-known 
(|la|. eq. 13.48 in Ref. [21]). They can be written in one 
formula as 



R 2 dz9 2 



dT 2 
1 



e - 2m/R 



dR 2 



(6.22) 



ds 2 = -dt 2 



R, r 



F 



dr 2 - 



-R 2 



dd 2 



1 • 2 

- sin 

e 



W~ri) dip 2 



(6.25) 

and the coordinates of (|6.25p are adapted to matching it 
to the planar or pseudospherical Szekeres metrics across a 
hypersurface of constant r. 

The matching requires that the intrinsic metric of a hy- 
persurface r = const and the second fundamental form of 
this hypersurface are the same for both 4-metrics. The 
match between the 3-metrics follows easily. Suppose the 
matching is done at r = b. The transformations to be ap- 
plied to (|6.25p are different for each value of e. With e = 
we transform the coordinates of (|6. 25|) as follows 



2\Jp' 2 + q' 2 , arctan(g'/p') 

S(b)(p-P(b),g-Q(b)) 
(p-P(b)) 2 + {q-Q(b)) 2 - 



(6.26) 



where to is a constant and e = ±1, 0. The metric with e = 1 
is the Schwarzschild solution; with e = and m > it is the 
Kasner solution in untypical coordinates, as is easy to ver- 



The transformation 
R = (9i 2 m/2) 1 / 3 
(2/(9m)) 1 f' i ^/x 2 +y 



to the well-known form 
T = (3/(4m)) 1 /3 Zi 

2 tp = tan -1 (x/y). 
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With e 

taah(i?/2) = 



1, we transform (|6.25p by 
1 



S(b) 
arctan 



q - Q(b) 



P - P(b) 



(6.27) 



After the transformation f|6 . 26|) . or, respectively, (|6.27p . the 
metric (I6.25P becomes 



ds 2 



-dt 2 
S(b) 



R. 



e + F 



dr 2 



Rf_ 



(dp 2 + dq 2 ) 



{p-P{b)) 2 + {q-Q(b)f 
(S(b)) 2 



+ e 



(6.28) 



In the single r = b hypersurface the 3-metric of (|6.28p has 
the same form as in (|2.ip . The two 3-metrics will coincide 
if their R(t,b) are the same at all times. This will be the 
case when 



M(b) = m , f(b) = F(b) , 



(6.29) 



since then both R-s obey the same differential equation, so 
it is enough to choose the same initial condition for both of 
them. The unit normal vector to the matching hypersurface, 
n a in the Szekeres metric is 

j^t RE' IE 

n Sa = (0, m, 0, 0) , m = — r —l , (6.30) 

V£ + / 

and in the vacuum metric (|6.28p it is n ea = 
(0, R' '/ 'y/e + /, 0, 0). In spite of these different forms, the 
terms (R' — RE' /E) and R' cancel out in the extrinsic 
curvature for each metric, and the only nonvanishing com- 
ponents of the second fundamental form of the r = b hy- 
persurface Kij are K22 = K33 = — V e + fEL/E 2 which are 
continuous across r — b by virtue of (|6.29p and (|6.28p . 

We see that, with the Szekeres mass function M being 
positive, the matching implies to > in both cases, and 
so the "exterior" vacuum solutions for e < are necessarily 
nonstatic. This, in turn, implies that any Szekeres dust 
model that matches on to them cannot be in a static state. 

Thus, in the most general case, the exterior metric for the 
planar Szekeres metric is the vacuum Kasner metric, and for 
the pseudospherical Szekeres metric it is the e = — 1 vacuum 
metric (|6.22p , both represented as in (|6.28p . 



VII. A FLAT MODEL OF THE SZEKERES SPACES 
WITH e / 

To visualise the geometric relations in the Szekeres spa- 
tial sections, we will construct an analogue of the Szekeres 
coordinate system in a flat 3-space. 



A. The quasi-spherical case. 



construct the appropriate coordinates in a plane. The setup 
will be axially symmetric, and after the construction is com- 
pleted we will add the third dimension by rotating the whole 
set around the symmetry axis. The foliating spheres inter- 
sect the plane along nonconcentric circles. The family of 
circles is such that their radii increase from to 00 while 
the positions of their centers move from the point (b, 0) to 
(+00, 0) in such a way that in the limit of infinite radius the 
circles tend to the vertical line x = 0. 

The family of circles is shown in the right half of Fig. |TTJJ 
it is given by the equation: 



V 



(7.1) 



where b is a constant that determines the position of the 
center of the limiting circle of zero radius which we will call 
the origin O, while u is the parameter of the family - the 
radius of the circles. 

For this family, we now construct a family of orthogonal 
curves. The tangents to the family (|7.ip have slope 



dy 
d.r 



y 



2 -x 2 +b 2 



2xy 



(7.2) 



which is the differential equation whose solution is (|7.1 
The orthogonal curves will obey the equation 



dy 
dx 



-2xy 



dx 



x 2 -y 2 - b 2 



(7.3) 



y 2 — x 2 + b 2 dy 2xy 

Note that this results from (|7.2j) by the substitution: 

(x,y,b) = (y',x',ib'). (7.4) 



Thus a solution of (|7.4|) results from (|7.ip by the same 
substitution and it is: 



+ (e\y — \/ v 2 — b 2 ^j = v 2 



(7.5) 



where e\ = ±1 and v is the parameter of the orthogonal 
family. As it happens, (|7.5|> is also a family of nonconcentric 
circles whose centres all lie on the y-axis, but the radius of 
the smallest circle is b. All the circles pass through the 
origin O at (x, y) — (b, 0). Fig. ITOl shows the x > part of 
both families. The double sign in (|7.5|) is needed to cover 
the whole right half of Fig. [T0J With only the + sign, only 
the y > sector would be covered. We did not include the 
corresponding double sign in (|7.ip because we wanted to 
cover only the x > half-plane with those circles. 

We now choose u and v as the coordinates on the plane 
and calculate the metric in these coordinates. From ()7.ip 
and (17.51) we find 



bv 



eie 2 - 



D ' 

D = v\J u 2 + b 2 + e2U\J v 2 



b 2 , 



(7.6) 



For the beginning, we will deal with the e = +1 case, 
i.e. with the foliation by nonconcentric spheres. We first 



where e 2 = ±1. The two solutions arise because, as seen 
from Fig. [TO] the pair of circles corresponding to a given 
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pair of values of (u, v) in general intersects in two points. 
The exceptional cases are u = (which corresponds to the 
single point (x, y) = (b, 0)) and v = ±b - when the v circle 
is mirror-symmetric in y, and the two intersection points 
have the same x coordinate. Using (j7.6j) we find 



dar + dy 2 



D 2 



b 2 



du 2 



v 2 — b 2 



dv 2 (7.7) 



To make the metric look more like Szekeres, we now 
transform the coordinate v as follows: 8 



b 2 + w 2 /A 



after which the metric (17.71) becomes 



(7.8) 



dx 2 + dy 2 



E = wD/b 2 



E 2 



(b 2 + w 2 /4y 

u 2 + b 2 



du 2 + u 2 dw 2 



\Ju 2 + b 2 - u+ (V« 2 + b 2 + uj 



w 

lb 2 ' 



(7.9) 



In the above, for a more explicit correspondence with 
the Szekeres solution, we have chosen e-i = —1 and 

\J (b 2 - w 2 /4) 2 = + (b 2 - w 2 /4), so that the term inde- 
pendent of w tends to zero as u —> oo. This will correspond 
to e — ► in the Szekeres metric. (The case e 2 = +1 results 
from (|7.9p by the inversion w = 4b 2 /w'.) 

By looking at the Szekeres metric (|2.ip we see that in 
(j7.9j) u simultaneously plays the role of r and of R. Let us 
follow the analogy. We are considering a flat 3-space (so 
far, only 2-plane). The 3-space t = const in the e = +1 
Szekeres metric will be flat when / = 0. Thus, if (j7.9j) is to 
become the metric of a flat space t = const in the e = +1 
Szekeres metric, then the coefficient of du 2 in (|7.9p should 
obey: 



(b 2 + w 2 /4) 2 
(u 2 + b 2 ) E 2 



-t ( E - uE. 

E 



(7.10) 



As can be verified, this holds. 

Now it remains to add the third dimension by rotating 
the whole configuration around the x axis of the initial 
coordinates. Thus, in {dx 2 + dy 2 ) we now treat y as a 
radial coordinate, we add <fi as the angle of the polar co- 
ordinates, and consider the metric (dx 2 + dy 2 + y 2 



We go back to (|7.6p and repeat the calculations with this 
3-dimensional metric. Thus, going to the Cartesian coordi- 
nates (y, z) = (y cos <j), y sin cf>) we thereby transform (w, (j>) 
to (p, q) — w(cos(j>, sin</>), or 



= vV +- 



4> = arctan(<7/p), 



(7.11) 



The transformation is a composition of two transformations: v - 
b/ sinx and \ = 2 arctan[(«i/(26)] . 




FIG. 10: A section through the family of nonconcentric spheres 
going over into a family of hyperboloids via a plane, that provides 
a model of flat space sections in Szekeres coordinates. The curves 
orthogonal to them are also shown. More explanation in the text. 



and after this the metric becomes 



d 53 2 



uE^ 



du 2 



E = Vi 



b 2 



\Ju 2 + b 2 + 



(p 2 + <7 2 ).(7.12) 



This is the axially symmetric (and flat) subcase of the 3- 
space t = const in the Szekeres metric of (|2.ip - (|2.5p 
corresponding to e = +1, B\ = B2 = P = Q = 0, S — 

2b 2 / (Vu 2 + b 2 + it) = 2 (Vi 



b 2 — u) and r — B = u. 



B. The Pseudospherical Case 

We can deal with the quasi-pseudospherical case in a sim- 
ilar way, but there is an important difference; in this case 
the 3-space of constant time can be flat only if its metric 
is pseudoeuclidean. This pseudoeuclidean space represents 
a space of Euclidean signature that has constant negative 
curvature. We will plot the coordinate curves in a Euclidean 
space, so it has to be remembered that they are distorted 
and do not represent the geometric relations faithfully. In 
particular, vectors or curves that are orthogonal in the pseu- 
doeuclidean metric will not look orthogonal in the plot. The 
space t = const can have a Euclidean signature, but then 
it must be curved. So we can say that we are representing 
the relations in a curved space with / 7^ by figures drawn 
on a flat plane, and hence the distortion. 

We begin with a family of right hyperbolae that fill the 
left half of the (x, y) plane and in the limit u — ► 00 tend to 
the straight line x = 0; at the end of the calculation we will 
rotate the whole collection around the x-axis. The family is 
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given by 



For the same reason as with (17.911, we have now chosen 



,x - ei 



\/u 2 — b 2 ) — y 2 —u 2 , 



(7.13) 



where t\ = ±1 9 and u > b is the parameter of the family. 
The family is shown in the left part of Fig. [TU] 

We again construct the family of curves orthogonal to 
these hyperbolae, but in the pseudoeuclidean sense. The 
hyperbolae (|7.13p solve the differential equation 



dy x 2 + y 2 + b 2 
dx 2xy 



(7.14) 



so the curves that are (pseudo) orthogonal to them obey 
the equation: 



2xy 



dy = 

da: x 2 + y 2 + b 2 ' 



(7.15) 



which is obtained from (j7. 14[) simply by interchanging x 
and y. We thus conclude that the solution is obtained from 
(|7.13p by the same interchange, and so it is 



— t2\/ v 2 — b 2 ) — x 2 = v 2 , 



(7.16) 



where e-2 = ±1. The two families (|7. 13|) and (|7. 15|) are 
shown together in the left half of Fig. [TO] 

By solving (|7. 13|> and ()7.16p for x and y we find: 



y = £263 



D] t = €\v\j u 2 — b 2 + e^uy/ v 2 — b 2 , 

where ei, e 2 £3 = ±1; from which we get 

, , o b 4 ( v 2 du 2 u 2 dv 2 
- dx 2 + dy 2 



i 2 — b 2 v 2 — b 2 



(7.17) 



(7.18) 



Again substituting for v with (|7.8|> . and adding the 
3rd dimension by rotating around the x axis in a 
similar way to (|7.1ip , the three dimensional metric 

(— dx 2 + dy 2 + y 2 dip 2 ") becomes 



b 2 



u 2 (dp 2 + dq 2 )} 



b 2 



\Ju 2 -b 2 + 
it 5 



(p 2 + q 2 ). (7.19) 



Just as for the spherical case, it can be verified that the 
Szekeres relation, analogous to (|7.10p . is obeyed. It reads 
here 



(u 2 ~b 2 )E 1 



Ei 



(7.20) 



9 This double sign is necessary in order that the hyperbolae fill the 
whole left half-plane. With only one sign they would fill only the 
part of the half-plane that lies to one side of the u = b curve. 



e 3 = -ei and ^(b 2 -w 2 /4) 2 = +(b 2 -w 2 /4). Then 
the sign of ei becomes irrelevant, since Dh oc ex, and only 
Dh appears in the metric. The metric (|7.20p corresponds 



to dZTp with e = -1, Bx = B 2 = P = Q = 0, S = 
2b 2 / (Vu 2 - b 2 + u)=2 (\/u 2 - b 2 - u) and r = R = u. 

Fig. [TUl shows the junction of the spaces of (|7.12p and 
(I7.19P - it represents the Szekeres t = const space con- 
sisting of nonconcentric spheres (right half of the picture) 
that tend to the plane x = from one side, and the family 
of hyperboloids (left half of the picture) that tend to the 
same plane from the other side. This shows how spheri- 
cal surfaces can go over into hyperboloidal surfaces within 
the same space. We repeat that only the right half of the 
picture faithfully represents the geometry of the flat space 
in coordinates defined by the spheres; the left half is a dis- 
torted image of either a curved 3-space or of a flat 3-space 
that has the pseudoeuclidean signature ( — h +)■ 

Note that the plane that separates the family of spheres 
from the family of hyperboloids has, in each family, the 
equation u — > 00 (to see this, solve (jZJj) and (|7.13p for x; 
in each case one of the solutions resulting when u — > 00 
is the plane x = 0). In this limit \Ju 2 ± b 2 — u — > 0, 
which, on comparing (|7.12p and (|7.19|) with (|2.2p , is seen 
to correspond to e = 0, just as it should. 



VIII. PHYSICAL DISCUSSION 



A. Role of R 



In the metric (|2.ip and in the area integral, A = 
R 2 J l/E 2 dpdq, the factor R 2 multiplies the unit sphere 
or pseudosphere, and therefore determines the magnitude 
of the curvature of the constant (t,r) surfaces. It is also 
a major factor in the curvature of the constant t 3-spaces. 
Therefore we view it as an "areal factor" or a "curvature 
scale". However, when e < it is not at all like a spherical 
radius. We note that when e = — 1, there can be no origin, 
but R can have maxima and minima as r varies, while in 
the e = case, R cannot have extrema, and it can only 
approach zero asymptotically. 



B. Role of M 

In (|2.5p , M looks like a mass in the gravitational poten- 
tial energy term of the evolution equation (|2.5p , while in 
(|2.6p M determines the deceleration of R. For e = +1, 
the function M(r) plays the role of the gravitational mass 
contained within a comoving "radius" r, but this interpre- 
tation is geometrically and physically correct only in the 
quasi-spherical model, where the surfaces of constant r are 
nonconcentric spheres enclosing a finite amount of matter. 
For e < however, R is not the spherical radius that is an 
important part of these ideas in their original form, and M 
is not a total (gravitational) mass, since the constant t & 
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r surfaces are not closed. Consequently these ideas need 
revising. 

The impossibility of an "origin" or locus where M and 
R go to zero when e = — 1 means that M must have a 
global minimum, and indeed regular maxima and minima 
in R and M are possible. Therefore the local M value is 
not independent of it's value elsewhere, and integrals of the 
density over a region always have a boundary term, sug- 
gesting the value of M (rather than its change between two 
shells) is more than can associated with any finite part of 
the mass distribution. In e = models, an asymptotic "ori- 
gin" is possible, but not required, and regular maxima and 
minima in R and M are also possible asymptotically. So, 
with an asymptotic origin (as occurs in the planar foliation 
of k = — 1 RW) the boundary term can be set to zero, but 
not with an asymptotic minimum in M and R. 

Nevertheless, the central roles of R and M are confirmed 
by the fact that the 3 types of Szekeres model can be joined 
smoothly to vacuum across a constant r surface at which 
the values of R and M must match (section |WB]). The vac- 
uum metric "generated" by the Szekeres dust distribution 
must have spherical, planar, or pseudospherical symmetry, 
and in each M is the sole parameter, while R is an areal 
factor. 

We note that, even in the Poisson equation, the grav- 
itational potential does not need to be associated with a 
particular body of matter, and indeed it is not uniquely de- 
fined for a given density distribution. 

Therefore we find that M is the mass factor in the grav- 
itational potential energy. 

C. Role of / 

As shown in section flTD] and as is apparent from the met- 
ric (|2.ip , the function / determines sign of the curvature 
of the 3-space t — const, as well as being a factor in its 
magnitude (c.f. [27]). In the quasi-spherical case, e = +1, 
this 3-space becomes Euclidean (represented in odd coordi- 
nates) when / = 0. In the quasi-pseudospherical case, with 
/ = it becomes flat but pseudoeuclidean: the signature is 
( — h +). In the quasi-planar case, the equations show the 
value / = is not possible, and thus the quasi-planar case 
does not in fact admit flat 3-dimensional subspaces. 

We also see from above that / appears in the gravita- 
tional energy equation as the total energy per unit mass of 
the matter particles, and we do not need to revise this in- 
terpretation. Therefore, this variable has the same role as 
in quasi-spherical and spherically symmetric models. 

D. Role of E 

We have seen in section [TH] that for e = +1, E'/E is the 
factor that determines the dipole nature of the constant r 
shells, and for e = —1, it is the pseudospherical equivalent 
of a dipole, except that the two sheets of the hyperboloid 
contain half the dipole each, and only one of them can 



be free of shell crossings. For e = 0, the effect of E'/E is 
merely to tilt adjacent shells relative to each other, but only 
the zero tilt case (E' = 0) is free of shell crossings. 

The shell separation (along the r lines) decreases mono- 
tonically as E'/E increases. If E' = it is uniform, other- 
wise it is minimum at some location and diverges outwards. 

For pseudospherical models, which must have / > 1, 
(|4~T4j> and ([2371) snow that if /7( 2 /) > M'/(3M) every- 
where and there are no shell crossings, the density is at all 
times monotonically decreasing with E'/E, but asymptoti- 
cally approaches a finite value as E'/E diverges. Therefore 
the density distribution on each shell is that of a void, but 
the void centres on successive shells can be at different (p, q) 
or (8, <fi) positions, in other words, the void has a snake-like 
or wiggly cylinder shape. The minimum density is only zero 
if M'/(3M) = -(E'/E) m3X . Far from the void, at large 
9, the density is asymptotically uniform with p &i q (i.e. 
with (/>), but can vary with r, though fairly gently compared 
with the void interior. If f'/(2f) < M'/(3M) everywhere 
, an initial void can evolve into an overdensity. Intuitively, 
it makes sense that there should be an initial underdensity, 
since too strong a tube-like overdensity would cause outer 
shells to expand much less rapidly, but this would cause shell 
crossings in models with hyperbolic evolution, / > 0. 

The location of the density minimum (or maximum) on 
each sheet is given by (|3.1ip and (|3.13f) . so their values 
are limited by the no shell crossings condition (j3. 15|) . Their 
rates of change depend on S" , P" and Q" and since the 
latter are not directly limited, they could be arbitrarily large 
at any one point, however (|3.15p implies that for any given 
r\ and ri 

2 f\p'P" + Q'Q")dr + {P[f + (Q[) 2 

<2 [ 2 S' S" dr + (S[) 2 (8.1) 
J n 

which means there is a limit on how far the location of the 
minimum can move for a finite change in r. The density is 
affected by E'/E at all times except near a simultaneous 
bang or crunch. 

IX. CONCLUSIONS 

We have analysed the Szekeres metrics with quasi- 
pseudospherical and quasi-planar spatial foliations, and es- 
tablished their regularity conditions and their physical prop- 
erties. 

For the quasi-pseudospherical case (e = —1), each con- 
stant r shell is a two-sheeted right hyperboloid (pseudo- 
sphere), each mapping to only part of thep-q plane, but only 
one sheet can be free of shell crossings, and only if E'/E 
has a negative maximum, going to — oo in the asymptotic 
regions of the sheet (where E — ► 0). The effect of E'/E 
can be called the pseudospherical equivalent of a dipole, but 
half the dipole is in the disallowed sheet. At this maximum 
the constant r shells are closest and the density has an ex- 
tremum — a minimum if (j4.14|) holds, otherwise it starts as 



20 



a minimum, but evolves into a maximum. Far from the the 
extremum, the density becomes uniform with p &i q, but can 
still vary with r. The location and value of the extremum 
depend on the derivatives of S(r), P{r) and Q(r), and so 
the density extremum can vary in magnitude and makes a 
wiggling, snake-like path. We also find that on a spatial 
section R can have extrema, but cannot be zero. The con- 
ditions for no shell crossing are weaker than for LT models, 
allowing R 1 or M' to become negative, though there is an 
extra condition relating the Szekeres functions S' , P' and 
Q' . (In contrast, for spherical foliations the no shell crossing 
conditions are not weaker than in LT.) 

We found the quasi-planar case (e = 0) was the hardest to 
understand. Only the plane symmetric case, E' = can be 
free of shell crossings, spatial sections can have neither zeros 
nor extrema of R, except asymptotically, and it isn't possible 
to make it spatially flat, / = (except in the Kantowski- 
Sachs-like limit [20]), so as a complete manifold it turns out 
to be the most restricted once physical regularity conditions 
are imposed. However a Szekeres spacetime can consist of a 
region with a quasi-spherical foliation joined across a planar 
boundary to a region with a quasi-pseudospherical foliation, 
as visualised in the the 3-d model of section [VTTI 

It was necessary to take particular care in analysing the 
meaning of R and M. Although the evolution of R obeys an 
energy equation with a term M/R that is very like a gravi- 
tational potential, R cannot be a spherical radius as the sur- 
faces it multiplies are not closed. Similarly, because there's 
always a boundary term when e = — 1, M is not solely 
determined by the matter inside a finite region, though its 
change in value between two constant r shells may possibly 
be associated with the matter bewteen them. 

Nevertheless, R is very closely tied to the curvature of 



the p-q 2-surfaces and to their areas, so it is a "curvature 
scale" or an "areal factor". The Poisson equation and the 
EFEs only relate field derivatives to the local matter, so in 
general the gravitational potential $ has no simple connec- 
tion to a volume integral of the density p. Similarly, in the 
planar and pseudospherical foliations studied here, we see 
an example in which the gravitational potential energy term 
in the R evolution equation (|2.5p is affected not only by the 
density and curvature in a finite region, but also by bound- 
ary values determined by the distant density and curvature 
distribution. We view M as a "potential mass" since it 
is a quantity with dimension mass that determines a grav- 
itational potential energy through M/R, and acceleration 
through M/R 2 . M is the key gravitational field parame- 
ter, and it relates the curvature scale R of the comoving 
surfaces to the potential energy of the R equation. 

Having understood some key features of the e < Szek- 
eres metrics, it should now be easier to construct useful 
models out of them. 



Acknowledgments 

The work of AK was partly supported by the Polish Min- 
istry of Science and Education grant no 1 P03B 075 29. AK 
expresses his gratitude for the Department of Mathematics 
and Applied Mathematics in Cape Town, where most of 
this work was done, for hospitality and perfect working con- 
ditions. CH thanks the South African National Research 
Foundation for a grant. This research was supported by an 
award from the Poland-South Africa Technical Agreement. 



[1] P. Szekeres (1975a) Comm. Math. Phys. 41, 55-64. 
[2] W. B. Bonnor, N. Tomimura (1976), Mon. Not. Roy. Astr. 
Soc. 175, 85. 

[3] S. W. Goode and J. Wainwright, Mon. Not. Roy. Astr. Soc. 
198, 83 (1982). 

[4] S. W. Goode and J. Wainwright, Phys. Rev. D26, 3315 
(1982). 

[5] P. Szekeres (1975b) Phys. Rev. D 12, 2941-8. 

[6] W.B. Bonnor (1976a) Nature 263, 301. 

[7] W.B. Bonnor (1976b) Comm. Math. Phys. 51, 191-9. 

[8] W. B. Bonnor, A. H. Sulaiman and N. Tomimura (1977), 

Gen. Rel. Grav. 8, 549-559. 
[9] M.M. de Souza (1985) Rev. Bras. Fiz. 15, 379. 
[10] W. B. Bonnor, Class. Q. Grav. 3, 495 (1986). 
[11] W. B. Bonnor, D. J. R. Pugh, South Afr. J. Phys. 10, 169 
(1987). 

[12] P. Szekeres, in: Gravitational radiation, collapsed objects 
and exact solutions. Edited by C. Edwards. Springer (Lec- 
ture Notes vol. 124), New York, p. 477. 

[13] K. Bolejko (2006), Phys. Rev. D 73, 123508. 

[14] C. Hellaby and A. Krasinski (2002) Phys. Rev. D 66, 
084011, 1-27. 

[15] J. Plebanski and A. Krasinski (2006) An Introduction to 



General Relativity and Cosmology, Cambridge U P, ISBN- 

13 978-0-521-85623-2. 
[16] A. Krasinski (1997) Inhomogeneous Cosmological Models, 

Cambridge U P, ISBN 521 48180 5. 
[17] G. F. R. Ellis (1967). Dynamics of pressure-free matter in 

general relativity, J. Math. Phys. 8, 1171. 
[18] M. Cahen and L. Defrise (1968) Comm. Math. Phys. 11, 

56. 

[19] W.B. Bonnor (1997) Mathematical Reviews 97 g 83 026 

(July 1997, p 4592). 
[20] C. Hellaby (1996) Class. Q. Grav. 13, 2537-46. 
[21] H. Stephani, D. Kramer, M. MacCallum, C. Hoenselaers 

and E. Herlt (2003) Exact Solutions of Einstein's Field 

Equations, Second Edition, Cambridge U P, ISBN 521 

46136 7. 

[22] B.K. Berger, D.M. Eardley and D. W. Olson (1977) Phys. 

Rev. D 16, 3086-9. 
[23] C. He llaby (1996) J. Math. Phys. 37, 2892-905. 
[24] http://www.maplesoft.com . 

[25] P. Musgrave, D. Pollney and K. Lake (2001), GRTensorll 
version 1.79, running under Maple V release 3. (Physics 
Department, Queen's University, Kingston, Ontario, K7L 
3N6, CANADA; E-mail: grtensor@astro.queensu.ca, web: 



21 



http://grtensor.org . ) 
[26] C. Hellaby and K. Lake (1985) Astrophys. J. 290, 381-7. 
[27] C. Hellaby (1987) Class. Q. Gray. 4, 635-50. 



APPENDIX A: THE E' = LOCUS AS A GEODESIC OF 
THE 2-D HYPERBOLOID 



The set E' = in the (p, q) surface of the metric (12. ip 
with e = — 1 is a geodesic in that surface. Proof: 

Calculate E' and rewrite the result in the {9,<f>) variables 

of 

S' cosh 8 + smh6(P' cos c/> + Q' sin <j>) 
~~ coshfl- 1 ' ' ' 

The solution of the equation E' = is 

5" 



tanh6* 



P' cos d> + Q' sin < 



(A2) 



Choose ^ as the parameter on the curve 6(4>) given by (|A2p . 
The tangent vector to this curve then has the components 

k a - (— i\ — - s'i-P'^^ + Q' 005 ^) 

~ \d(f>' J ' d0 ~ (P'cos0 + Q'sin0) 2 - S" 2 ' 

(A3) 

The metric (|2.38p , and its nonzero Christoffel symbols, in 
the (9,<p) coordinates, are 



ds 2 2 = d9 2 + smh 2 6>d<^ 2 , 

1 1 = -sinh0cosh<9, I z f.=coth6>. (A4) 

22 1 I 1 2: I 



Thus the equations of a geodesic are 

d 2 6> d(9 

— -77 — sinh 6* cosh 9 = A—, 
d0 2 a<p 

d.9 

2coth6»— = A, (A5) 



where A is an unknown proportionality factor. The second 
equation above defines A, which is 

2(-P'sin0 + Q'cos</>)(P'cos0 + Q'sin0) , Ac . 

A = n 1 (A6) 

(P'cos0 + Q'sin0) 2 - S' 2 

and then the first of (|A5P is easily verified using (|A3P and 

tanh 9 



S'{P' cos0 + Q'sin 0) 
(P'cos(/>-t-Q'shi0) 2 - S' 2 



(A7) 



□ . 



APPENDIX B: THE BIANCHI TYPE OF THE 
PSEUDOSPHERICAL VACUUM MODEL. 



sinh 9 cosh 9 = 



1 - tanli 9 



Lower indices will label vectors, the upper indices will 
label the coordinate components of vectors. The Killing 
vector fields for the metric f|6.22|> with e = — 1 are: 

c a i a c a 

1 =01, K 4 = 3 , 

k^* = costp82 a — cothf9shi<^o~3 Q , 

k 3 a = sin^ 2 Q +coth?9cos</J(53 Q! . (Bl) 

The commutators are: 

[fci.fci] = 0, 1=1,2,3, 

[fe, fc 3 ] = ki, [fc 2 , fc 4 ] = fc 3 , 

[*3,*4] = -fe- (B2) 

The Bianchi algebra must thus include k\ and a 2- 
dimensional subspace of {k 2 , fc 3 , £4}. Consequently, out of 
the set {^2,^3,^4} we choose two linear combinations, I 
and m, that span a 2-dimensional Lie algebra, i.e. have the 
property [£, m] = a£ + f3m. This can be done in many ways; 
one example of such a combination is 

£ = k 2 , m = fc 2 + fc 3 + fc 4 , (B3) 

for which we have 

[£, m]=m-£. (B4) 

This is not a standard Bianchi basis. To obtain a standard 
basis (see Ref. [15]) we take such combinations of k\, I 
and m that are equivalent to 

w\ = 2k 2 + k 3 + fe 4 , w 2 — k\ — fc 3 — fc 4 , 

W3 — k\ + fc 3 + fc 4 . (B5) 

The commutation relations are now 

[wi,w 2 ]=W2-w 3 = [ws,Wi], [w 2 ,w 3 ]=0, (B6) 

and this is the standard form of the Bianchi type III algebra. 



